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ABSTRACT 

A detailed analysis is presented of 33 RR Lyrae stars in Pisces observed with the 
Kepler space telescope over the 8.9-day long K2 Two-Wheel Concept Engineering 
Test. The sample includes not only fundamental-mode and first overtone (RRab and 
RRc) stars but the first two double-mode (RRd) stars that Kepler detected and the 
only modulated first-overtone star ever observed from space so far. The precision of the 
extracted K2 light curves made it possible to detect low-amplitude additional modes 
in all subtypes. All RRd and non-modulated RRc stars show the additional mode at 
Px/Pi ~ 0.61 that was detected in previous space-based photometric measurements. 
A periodicity longer than the fundamental mode was tentatively identified in one 
RRab star that might belong to a gravity mode. We determined the photometric 
[Fe/H] values for all fundamental-mode stars and provide the preliminary results of our 
efforts to fit the double-mode stars with non-linear hydrodynamic pulsation models. 
The results from this short test run indicate that the K2 mission will be, and has 
started to be, an ideal tool to expand our knowledge about RR Lyrae stars. As a 
by-product of the target search and analysis, we identified 165 bona-fide double-mode 
RR Lyrae stars from the Catalina Sky Survey observations throughout the sky, 130 
of which are new discoveries. 
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1 INTRODUCTION 

RR Lyrae stars are pulsating low-mass horizontal-branch 
stars. Since they can be found at the intersection of the 
horizontal branch and the classical instability strip, a rel¬ 
atively small region of the Hertzsprung-Russell-diagram, 
their absolute brightnesses are fairly similar, making them 
good distance indicators. The large number of RR Lyrae 
stars make them good tracers of h alo structures around 
the Milky Way (iSesar et all l2013l l and even in nearby 
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galaxies (e.g. iBrown et al.l 120041 : ISaraiedini et al.l 1 20061 : 
lHaschke. Grebel fe DuffatJ 20121 1. The period ratios of the 


rare double-mode stars can be used to determine the masses 
of RR Lyrae stars independent from stellar evolution theory. 


The Kepler space telescope observed ~ 45 RR Lyrae 
stars in the Lyra-Cygnus field of the prime mission, spanning 
a wide metallicity and brightness range, from the eponym of 


the class, RR Lyr (Kp = 7.862 mag. lKolenberg _et_alj 


down to about Kp = 17.5 magnitudes Nemec et al.1 


2011 ), 


201311 . 


Although the Kepler mission was designed t o detect su btle 
light changes caused by planetary transits dBorucki et al.1 
I 2 OIOII and the high-amplitude variations of the RR Lyrae 
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stars required additional processing, space-based photome¬ 
try proved to b e essent i al to furth er our understanding of 
RR Lyrae stars f Benko et alJl201 4'). 

Kepler observed mostly fundamental-mode stars, with 
only four first-overtone and no double-mode stars in the 
sample. Therefore most of its findings have been limited 
to the RRab cl ass. The results include the discovery o f 
period doubling dKolenberg et al.l 120101 : ISzabo et al. 2010|b 
a we a lth of low-amplitude addit ional modes ( Benko et al.1 
120101 : iGuggenberger et al.1 l2012l l and the discovery that 
most Blazh k o star s show multiple or irregular modulations 
llBenko et abl^OM f. The observations also revealed that the 
new phenomena (period doubling and a dditional mod es) ar e 
strictly limited to the modulated stars dNemec et al.ll201l|j . 
and a recent analysis of stars observed by th e C 0 R 0 T s pace 
telescope confirmed this dichotomy dSzabo et al.ll2014h . In 
contrast, all four first-overtone stars observed by Kepler 
(and two more in the C 0 R 0 T sample) exhibit the same 
additional mode at the frequency ratio / 1//2 ~ 0.60 — 
0.64, without any signs of modulation in the main mode 
dMoskalik et al.1 l2013il . 120151 ') . The origin of this additional 


mode is uncertain. We note that we refer to this mode as 
fx instead of /2 to avoid any confusion with the suspected 
second-overtone signals in the RRab stars. 

The primary mission of Kepler gathered four years of 
(quasi)-continuous observations. These allowed for more spe¬ 
cific studies, such as following the very unexpected vanish ing 
of the Blazhko effect in RR Lyr (|Le Borgne et al.l I2014T) or 
the first att empt to search for c haos in the modulation of an 
RRab star dPlachv et al.ll20141 ). and ev en the search for ar - 
tificial origins in the period variations dHippke et al.ll2014h . 
After the unfortunate failure of the second reaction wheel, 
however, maintaining the original pointing was not feasible 
any more. 

The new I<2 mission of the Kepler space telescope uses 
orientations along the Ecliptic where the space craft can bal¬ 
ance against the radiation pressure of the Sun dHowell et al.l 
120141 '). This new strategy results in short, ~ 75 day long ob¬ 
serving campaigns that scan various fields in the sky, reach¬ 
ing targets the original field lacked. In this paper we detail 
the results from the very first measurements of the K2 mis¬ 
sion. Section 2 describes the K2-E2 run and the extraction 
of the light curves. Results are detailed in Sections 3 and 4 
for RRd-RRc and RRab stars, respectively. Conclusions and 
future directions are summarized in Section 5. 


2 K2 DATA AND REDUCTION 

Kepler observed a stellar field around the vernal equinox 
point in Pisces (center coordinates: a = 359°, 8 = —2°) 
between 4th-13th February 2014. The primary goal of this 
K2 Two-Wheel Concept Engineering Test (hereafter K2-E2) 
was to test the performance of the telescope in fine guidance 
mode. As well, the observations of nearly 2000 targets were 
made available for the scientific community. Instead of the 
tight target pixel masks used in the primary mission, larger 
subframes were stored around the targets in order to avoid 
any flux loss due to the drift of the telescope. A ~ 1 pixel 
variation in the position of the stars was observed during the 
8.9-day long K2-E2 data. The boresight was shifted by mul¬ 
tiple pixels 2.3 days into the test, at BJD = 56695.36, when 


Table 1. Sample of the complete K2-E2 light curve file. The full 
table is accessible online. 


EPIC 
type: RRd 

BJD-2450000 

(d) 

Kp 

(mag) 

A Kp 
(mag) 

60018653 

1860.0503350 

13.90840 

0.00054 

60018653 

1860.0707686 

13.92217 

0.00054 

60018653 

1860.0912023 

13.92438 

0.00054 

60018653 

1860.1116358 

13.91703 

0.00054 

60018653 

1860.1320694 

13.89827 

0.00053 


the telescope locked successfully into fine guidance mode. 
Fortunately the 50x50-pixel target masks were more than 
enough to contain these movements. 

We identified 34 potential RR Lyrae stars in the K2- 
E2 sample and extracted their photometric data with the 
PyKE softwar^], developed for the Kepler mission b y the 
Kepler Guest Observer Office dStill fc Barclay! l2012| j . All 
targets were observed in long cadence mode, with a sam¬ 
pling of 29.4 minutes. To counter the effect of drifting and 
capture all the flux, we defined relatively large pixel aper¬ 
tures for the stars, with 1-2-pixel-wide halos around the ac¬ 
tual point spread functions. We applied two different masks 
for the positions before and after the boreshight shift. We 
used the background-corrected SAP (Simple Aperture Pho¬ 
tometry) fluxes calculated by the kepextract routine. The 
photometric precision spreads between 3.5 /rmag to 7 mmag 
per data point for a 10.5 and a 17.4 mag star, respec¬ 
tively. In some cases we had to apply small amounts of 
shift and/or scaling to the data to connect the two parts 
of the light curves. These differences were likely caused by 
different pixel sensitivities in the two positions. Some stars 
are affect ed by t he Fin e G uidance Sensor clock crosstalk 
patterns (ivan Cleeve fc Caldwellll2009li . The faintest star in 
the sample, EPIC 60018780, suffers from extensive negative 
video crosstalk from a bright target elsewhere on the same 
module: dark CCD columns cross the image of the star. As 
the image of the star moves in and out of the dark column, 
jumps up to 0.1-0.2 mag appear in the light curve of this 
star, especially in the first 2.3 days. 

Without the third reaction wheel, radiation pressure 
from the Sun slowly turns Kepler away from the desired 
spacecraft attitude, and rolls it about t he opti cal axis 
(the X-axis of the spacecraft, iHowell et al] l2014|j . These 
disturbances are periodically corrected with the on-board 
thrusters. For some stars the frequency of the manoeuvres, 
/corr = 8.16 d” 1 , and its harmonics can be detected in the 
Fourier spectra. (We note that the K2-E2 run tested a 3- 
hour adjustment frequency, but later campaigns use 6-hour 
intervals, corresponding to f co rr = 4.08 d _1 .) Rolling of the 
field-of-view means that stars close to the center are least 
affected but targets close to the edges can experience con¬ 
siderable drifting. Although these peaks are well-defined and 
can be easily excluded from the analysis, they indicate that 
special care is needed to define the target apertures and to 
correct for pixel sensitivity variations. 

Undoubtedly, there will be better and more precise 
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Table 2. Sample of the frequency tables of the entire K2-E2 data set: frequencies, amplitudes, phases, the respective uncertainties 
and the identifications of the peaks (/o: fundamental mode, f\: first overtone, fx , / 2 , fg- additional modes). Super-Nyquist frequency 
components are marked with asterisks. The entire table is accessible in the online version of the paper. 


Star 

F (d- 1 ) 

A (Kp mag) 

<p (rad/27r) 

A F (d" 1 ) 

A A (Kp mag) 

A (f) (rad/27r ) 

ID 

60018653 

1.85382 

0.04641 

0.92372 

0.00031 

0.00040 

0.00080 

/o 

60018653 

2.48564 

0.13988 

0.82381 

0.00010 

0.00035 

0.00026 

h 

60018653 

3.7076 

0.00422 

0.2164 

0.0034 

0.00034 

0.0088 

2/o 

60018653 

5.561 

0.00105 

0.573 

0.014 

0.00037 

0.035 

3/o 

60018653 

4.97128 

0.02273 

0.1761 

0.00063 

0.00036 

0.0016 

2/r 


ways to extract K2 photometry as the mission evolves. 
However, (semi-) automated methods are usually optimised 
for relatively quiet stars and not for large-amplitude, 
sh ort-perio d pulsators. One such pipeline, developed by 
IVanderburg fe Johnsonl (120141 ') . for example, admittedly 
failed to improve the light curves of RR Lyrae stars and 
eclipsing binaries. Therefore we were satisfied with the level 
of accuracy we achieved with the above method for this ini¬ 
tial study. The light curves (along with the frequency ta¬ 
bles, see Section [2.11) are available online for further analy¬ 
sis alongside the paper and at the Kepler Investigations at 
Konkoly webpagqj. Samples are provided in Tables [Q and [5] 

2.1 The K2-E2 RR Lyrae sample 

Prior to the K2-E2 test, a call was issued within KASC (the 
Kepler Asteroseismic Science Consortium) to provide tar¬ 
gets within the approximate position of the field in Pisces. 
We proposed 143 RR Lyrae stars, largely based on the R Rab 
catal ogue of the Catalina Sky Survey (CSS, iDrake et all 
l2013al a) but also on the ASAS database and the SIMBAD 
identifications. By the time the K2-E2 data became avail¬ 
able, the detaile d classifications of periodic variables in CSS 
were published dDrake et al.ll2014h . allowing for a follow-up 
search in the target masks. We found 34 stars at or close 
to the positions of the 1952 K2-E2 targets based on these 
catalogues. Comparison with the proposed list revealed that 
several of the originally proposed targets ended up off sili¬ 
con, while others were not selected for observation. 

One of the stars, EPIC 60018241 (A SAS J234439 - 
0148.6) turned out to be an eclipsing binary (I Conroy et al.l 
120141) . but the other 33 stars are RR Lyrae stars. Surpris¬ 
ingly, we identified RR Lyrae stars from all subclasses, in¬ 
cluding two classical double-mode (RRd) stars and a po¬ 
tential strongly modulated first-overtone (RRc) candidate. 
Basic parameters, i.e. the EPIC {Kepler Ecliptic Plane In¬ 
put Catalog) identification number, coordinates, brightness, 
period, other identifications, and subclass, are summarised 
in Table [3] for RRd and RRc stars. The parameters of RRab 
stars are listed in [5] 

We use d the Period04 softwa re to calculate the Fourier 
transforms llLenz fe Bregerll2005l) . The frequency tables of 
the K2-E2 sample are summarised in Tabic [5] Here we in¬ 
cluded only the resolved frequency peaks, i.e. unresolved 
residuals that indicate amplitude and/or phase variation are 
not listed. In some cases high-order harmonics reflect back 

2 http://konkoly.hu/KIK/data.html 


from the Nyquist frequency (/jv = 24.4684 d _1 ), but we fit¬ 
ted the super-Nyquist values in all cases. These are marked 
with asterisks in the frequency table. 


3 DOUBLE-MODE AND FIRST-OVERTONE 

STARS 

3.1 Double-mode stars 

We identified two RRd-type stars in the field. Both stars, 
EPIC 60018653 (CSS J235643.2+030119) and 60018662 
(CSS J234207.3+003253), were misclassified in the Catalina 
sample, as RRc and RRab stars, respectively, but the K2 
light curves revealed that they are actually classical double¬ 
mode RR Lyraes (Figure [T]). These are the first RRd stars 
observed by Kepler and only the third and fourth from 
space. AQ Leo, the brightest of th e class, was observed 
by MOST (iGruberbauer et al.ll2007h . and the CoRoT data 
of another star, CoRoT ID 0101368812, was analysed by 
IChadidl (I2012T ) and ISzabo et al.l (120141) . 

3.1.1 Period ratios and physical properties 

The two K2 stars have similar pulsation properties. Their 
fundamental-mode periods are Po = 0.53943 d and Po = 
0.55900 d, while the period ratios are Pi/Po = 0.74581 and 
Pi/Po = 0.74678 for EPIC 60018653 and 60018662, respec¬ 
tively. These values put them into the long-period end of the 
Petersen diagram that displays the period ratios against the 
fundamental-mode period of the stars (Figure [2]). In fact, 
EPIC 60018662 is near to the long-period edge of the area 
covered by double-mode stars. 

We constructed a Petersen diagram to put the two 
stars into broader context. The diagram in Figure [5] in¬ 
cludes the largest sample of field RRd stars constructed 
so far. One of us (RP) searched for bona-fide double-mode 
stars i n the Cata lin a Su rveys Catalog of Periodic Variable 
Stars (IDrake et al.l 120141) and identified 130 new RRd vari¬ 
ables. The details of this search and the list of stars are 
presented in Appendix m We also included the stars iden- 


tided in the LINEAR (Poleski 2014 

), ASAS/NSVS (IWils 

20051; Wils, Lloyd & Bernhard 2006[" 

Bernard & Wils 

2005 

Szczvgiel & Fabrvckv 20071), and Super WASP (Wils 

20 im 


databases. Stars from the Sagittarius dwarf galaxy and the 
OGLE bulge su r vey provide the backg round (grey crosses; 
ICseresniea l200ll ; ISoszvnski et al.l 120141 ') . It is interesting to 
note that even with the CSS sample, field stars outside the 
bulge seem to be restricted to Po ^ 0.45 d and Pi/Po ^ 
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Figure 1. Light curves (top row) and Fourier-spectra of the two RRd stars in the sample. Middle row: the original Fourier spectra. 
Fundamental-mode and first-overtone frequency peaks are labelled with fo and /i. Bottom row: residual spectra after prewithening the 
data with fo, /i and their linear combination peaks. The additional mode f \ and its significant subharmonic peaks and combinations 
are labelled. We also indicated the marginally detected 1/2 fx peak in EPIC 60018662. 


Table 3. Double-mode and first-overtone RR Lyrae stars observed during the Two-Wheel Concept Engineering Test run of the K2 
mission. Ol and FM indicate the periods of the first overtone and fundamental mode, respectively. We provide estimated Kp magnitudes 
for two stars (in italics). The brightness of EPIC 60018238 is too high in the EPIC (Kp = 9.83), compared to the flux values. The last 
star has no EPIC ID but was observed within the pixel mask of the engineering target EPIC 60042292. ^Modulation period. 


EPIC ID 

RA (deg) 

Dec (deg) 

Kp mag 

Ol period (d) 

Name 

FM period (d) 

CSS ID number 

60018653 

359.18015 

+3.022110 

13.760 

0.402311 

CSS J235643.2+030119 

0.539427 

1104128004573 

60018662 

355.53058 

+0.547972 

14.216 

0.417448 

CSS J234207.3+003253 

0.558996 

1101127011298 

60018224 

358.53385 

+0.963635 

10.600 

0.306156 

EV Psc 

— 

— 

60018238 

354.87500 

-3.631667 

12419 

0.274804 

ASAS J233930—0337.9 

- 

- 

60018678 

352.16375 

+1.028500 

14.862 

0.432426 

CSS J232839.2+010142 

- 

1101126021018 

(60042292) 

359.42500 

-1.839444 

14.502 

0.300110 

CSS J235742.1-015022 

17.07** 

1001128017680 


0.742, whereas the Sgr and bulge population extends dow n 
to Po ~ 0.35 d and Pi/Po ~ 0.726 llSoszynski et alJEoill) . 

The precise period values of pulsation modes in a star 
depend strongly on the fundamental parameters: mass, lu¬ 
minosity, effective temperature, and metallicity of the star. 
Therefore the simultaneous detection of two radial modes 
in one target can put strong constraints on its physical pa¬ 
rameters. This is especially important for mass determina¬ 
tion, as we currently have no other means to measure RR 
Lyrae masses. Although a f ew promisi ng binary candidates 
were recently discovered bv lHaidu et al.1 (120151/ . dynamical 
masses will take several years if not decad e s to m easure. Ac¬ 
cording to the calculations of lSzabo et al.l (|2004l ). RRd stars 
are confined to a narrow parameter range. We carried out a 
pilot study to see whether we can find good hydrodynamic 
models effectively. 

First we estimated the effective temperature, mass, lu¬ 
minosity and metal content ranges based on the pulsation 


model results of lSzabo et al.l (120041) . Compared to their Fig¬ 
ure 6/a, the period parameters of the two stars fall well be¬ 
low the Z = 0.001 models, therefore we selected two OPAL 
opacity tables with Z = 0.0001 and 0.0003 metal contents. 

Then a grid of non-linear hydrodynamic models was 
computed with the Florida-Budapest t urbu l ent convective 
code ( detail s of the code ar e prese nted bv lKollath fc Buchlerl 
(l200lJ) and iKollath et al.l (120021 )). Luminosity was varied 
between L = 48 — 60 Lq, and two temperature and mass 
ranges were selected: one between T e // = 6480-6540 K, 
M = 0.72-0.80 Mq, and the other between T e // = 6410- 
6470 K , M — 0.76-0.84 Mg. Increments were 2 or 3 Lq, 10 
K and 0.2 Mg. 

Models were iterated for 2000 fundamental-mode cycles 
and then the periods of both modes were calculated. We 
hand-picked 8 models close to the periods of EPIC 60018653 
and 6 models close to 60018662 and iterated them for a 
further 2000 cycles. After that, 7 models converged slowly 
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Figure 2. Period ratio versus the period of the fundamental mode 
of the RRd stars, also known as the Petersen diagram. Small 
grey crosses: OGLE bulge and Sgr dwarf stars, their sequence 
continues to even smaller values. Black squares: the two K2 stars; 
orange dots: stars detected in the LINEAR, ASAS, NSVS and 
SuperWASP surveys; blue circles: our CSS sample from Appendix 
S grey squares: CoRoT ID 0101368812 (left) and AQ Leo (right). 
Multiple identifications were not filtered. References for the data 
sources can be found in Section 13.1.11 

to single-mode pulsation and only 7 double-mode models 
remained. Six of those models are actually pairs that differ 
only in th e temperature step . The results agree with the 
findings of lSzabo et al.l (l2004ll that non-linear RRd models 
exist only in narrow parameter ranges. 

As Figure [3] illustrates we could not fit both stars with 
the same metal content value. For EPIC 60018653, two mod¬ 
els at T e ff = 6410-20 K, M = 0.76 Mq, L = 48 Lq and 
Z = 0.0003 are relatively close to the star on the Petersen 
diagram. Two models fall very close to EPIC 60018662 in 
this parameter plane: T e // = 6430-40 K, M = 0.78 Mq, 
L = 52 Lq and Z = 0.0001. We included, for comparison 
purposes, three more models that have similar fundamental¬ 
mode periods but different period ratios than the stars. De¬ 
tails of all seven models are included in Table [4] 

The mass values appear to be quite high: pulsation 
masses for RRab stars in th e origina l K eple r field, for exam¬ 
ple, fell below 0.66 Mq (iNemec et al.1 [201ll l. Unfortunately, 
in the absence of proper dynamic masses, all mass estimates 
are model-dependent. Low-mass linear model calculations 
can be fitted to the w hole range of the Petersen diagram 
l|Soszvnski et al.ll201ll ). However, non-linear RRd models in 
the Florida-Budapest code appear only at higher masses, 
typically above 0.7 Mq and they are still consistent with 
evolutionary models. We also note that the ability of current 
non-linear pulsation codes to properly model double-mode 
pulsat ion is still a matter of de b ate; s ee, e.g. iKollath et al.l 
d2002ll and ISmolec fc Moskali^ (120081 ) for opposing views. 
We also note that the lack of the fx mode in the mod¬ 
els does not invalidate our findings. If we assume that it is a 
non-radial mode, it could naturally arise in higher-dimension 
models with same the fundamental parameters. 

We emphasize that these are just preliminary results 
from a pilot study. Calculating double-mode non-linear mod¬ 
els is time-consuming, mainly because mode amplitudes may 
change very slowly. Nevertheless, we plan to carry out a more 


Table 4. Hydrodynamic model results for the two RRd stars, 
numbers correspond to EPIC 60018653 and 60018662. We also 
included the models where the first-overtone periods do not fit 
the observations, as indicated in Figure [3] 


Star 

Mass 

(Mq) 

Lum. 

( L o) 

Te f f 
(K) 

Z 

Po 

(d) 

P 1 /P 0 



Good model fits 



653 

0.76 

48 

6410 

0.0003 

0.54194 

0.74538 

653 

0.76 

48 

6420 

0.0003 

0.53906 

0.74545 

662 

0.78 

52 

6430 

0.0001 

0.56036 

0.74665 

662 

0.78 

52 

6440 

0.0001 

0.55737 

0.74674 


Models with different period ratios 


(653) 

0.76 

48 

6410 

0.0001 

0.53917 

0.74689 

(662) 

0.78 

51 

6410 

0.0003 

0.56031 

0.74516 

(662) 

0.78 

51 

6420 

0.0003 

0.55731 

0.74524 


0.7475 

0.747 


Nk: 0 

° 0 60018662 


0.7465 


. 

0.746 

_ 




0 

0 ° ° g o60018653 0 ° Q 


0.7455 

0 

° • # 



0 

* • % „ 


0.745 

0 0 0 



0.7445 

0 


- 

0.744 

0 0 




0.5 0.51 0.52 0.53 0.54 0.55 0.56 0.57 0.58 

P 0 M 

Figure 3. Hydrodynamic model results compared to the two RRd 
stars (blue squares) on the Petersen diagram. The small circles 
are the CSS RRd stars from Appendix [A] Large red dots are the 
Z = 0.0003 models, large crosses are the Z = 0.0001 models. 

extensive survey after more RRd light curves will become 
available from the K2 mission. 

3.1.2 Additional modes in RRd stars 

The stars AQ Leo and CoRoT ID 0101368812 both showed 
additional modes that could not be explained by linear com¬ 
binations of the fundamental mode and the first overtone. 
Comparison with the Fourier-parameters we derived from 
the K2 stars revealed striking similarities to these double¬ 
mode stars. One particular mode that falls to the frequency 
ratio of fi/fx ~ 0.615, was identified in both stars, and 
it has been recently identified i n a modulat ed RRd star 
from the OGLE database (ISmolec et al.l [20151) as well. Fig¬ 
ure [l] shows that we detected the same additional mode 
at frequency ratios fi/fx = 0.6162 in EPIC 60018653 and 
fi/fx = 0.6166 in 60018662. The same mode was reported 
not only in RRd stars but in RRc stars as well (Section 13. 21) . 

We detected the subharmonics, or half-integer peaks 
of the fx mode (1/2 and 3/2 fx) in EPIC 60018653. The 
presence of these frequencies indicates that the mode ex¬ 
periences period doubling. The peaks are broadened and 
their centres are slightly off from the exact half-integer val¬ 
ues at 0.499 fx and 1.495 fx, respectively. Temporal vari- 
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Figure 4. Folded light curves and Fourier-spectra of the three non-modulated RRc stars. Original spectra are shown with thin grey lines; 
the black thick lines are the residuals after prewithening with the nf\ (n = 1, 2,... 5) frequencies. The fx peak and the corresponding 
half-integer peaks are indicated. / corr is the instrumental signal of the correction manoeuvres. Labels of low-significance peaks (S/N ratio 
between 3 and 4) are in parentheses. 


ations in period doubling created forests of peaks in the 
spectra of RRab stars observed in the nominal mission and 
the highest-amplitude peaks w ere also displaced from the 
exact values (ISzabo et alj|2010h . Temporal variations have 
been detected in the fx mode and its sub harmonics in the 
Kepler RRc sample (IMoskalik et al.ll2015h . It is very likely 
that the broadened peaks in the K2-E2 sample are caused by 
temporal variations in the period doubling of the fx mode 
as well. The 1/2 fx peak is marginally detectable in EPIC 
60018662 too, but only with a signal-to-noise ratio of 2.5. 

iGruberbauer et~ahl (|2007l i observed the same 1/2 fx 
peak in AQ Leo, but identified this peak tentatively as the 
parent freq uency, an d the true fx additional mode as 
2 fi = fa. IChadidl J2012I) reported two other additional 
modes (/ 3 , fi) in the CoRoT star. However, the frequency 
of their fi — fi combination peak is equal to 1.49/x. If we 
accept this peak as a subharmonic, then fi itself becomes 
a combination peak (1/2 fx — fi) instead of an indepen¬ 
dent mode, making the star more similar to the rest of the 
RRc sample. Similar combinations can be detected in EPIC 
60018653 too. We can conclude that not only the fx addi¬ 
tional mode, but its period doubling effect is observable in 
all four RRd stars measured by space photometric missions 
so far. 


3.2 First-overtone stars 

Four RRc stars were identified in the K2-E2 data, the same 
number as observed in the original Kepler held. Three of 
those stars are normal first-overtone pulsators. The fourth 
star, CSS ,1235742.1-015022 is very different from the others: 
it exhibits strong amplitude and phase variations. 

The three non-modulated stars have different pulsation 
periods and light-curve shapes, as illustrated in Figure [4] 
EPIC 60018238 has the shortest period (P = 0.27474 d) 


and a very pronounced hump feature just before maximum 
light. We detected two different additional modes with very 
small amplitudes. The star seems to exhibit the ubiquitous 
fx mode, although we detected two, closely-spaced peaks 
at fi/fx = 0.6154 and fi/f' x = 0.6048 instead of a single 
frequency, with low significances (S/N ratios are ~ 3.3). We 
also identified the drift-correction frequency at 8.16 d _1 in 
the star. 


The light curves of the other two RRc stars, EV Psc 
(EPIC 60018224) and EPIC 60018678 show almost no shock- 
related features (humps or bumps). The pulsation period 
of EV Psc exhibits strong beating with the 1766-second 
long-cadence sampling period, so it is not covered contin¬ 
uously (Figure [4]). We detected the fx mode in both stars 
at frequency ratios of fi/fx = 0.6146 (EV PSc) and 0.620 
(60018678). The additional mode is very pronounced in EV 
Psc: the scatter of the folded light curve in Figure[4]is caused 
primarily by the mode itself. Moreover, we also detected the 
half-integer frequencies (1/2 and 3/2 fx) in EV Psc, corre¬ 
sponding to period doubling in the fx mode. Half-integer 
peaks can be located but are marginal in EPIC 60018678. 
It seems that the (period-doubled) fx mode is intimately 
connected to the first radial overtone: it was observed not 
only in the K2 RRc stars but in the original Kepler sample, 
in tw o CoRoT stars, and t he RRd stars desc ribed above as 
well (IMoskalik et al.ll2015l : ISzabo et al.ll2014l '>. 


In fact the fx mode has been detected in first-overtone 
Ceph eids Erst, from the data collected by the OGLE sur - 
veys (IMoskalik fc Kolaczkowskill2008l : IMoskalik et alJl2013h . 
More recently, the fx mode has been identified in 147 RRc 
stars in the OGLE bulge sample and in 14 RRc and 4 
RRd stars in the globular clus ter Messier 3 , respectively, 
further conhrming its ub iquity llNetzel. Smolec fc Moskalikl 
120151 : lJurcsik et al.ll20f.5ll. Y et the origin of the mode is still 
unclear: Dziembo wski ( 20120 hypothesised that high-degree 


































































33 RR Lyrae stars observed with K2-E2 7 


EPIC 60018224 


EPIC 60018238 


EPIC 60018678 


0.76 

— 0.72 
x 
-©- 

0.68 



4 8 12 16 20 

A x (t) (mmag) 


0 2 4 6 8 

A x (t) (mmag) 


0 2 4 6 8 10 
A x (t) (mmag) 


Figure 5. Time variation of the amplitude and phase of the fx 
in the non-modulated RRc stars (also known as loop diagrams). 
Arrows denote the direction of progression. Variation can be de¬ 
tected in all three stars: EPIC 60018224 (EV Psc) shows partic¬ 
ularly strong decrease in the amplitude of fx • 


/-modes (JL = 42 — 50) could explain the period ratios but 
those would require quite high intrinsic mode amplitudes. 
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3.2.1 Variation of the additional modes 

Beside the double peak in EPIC 60018238, we observed 
residual power in the Fourier spectra of both EPIC 60018224 
and 60018678 after prewithening with fx - These findings in¬ 
dicate that the fx mode undergoes a mplitude and/o r p hase 
vari ations, similar to t he findings of Mosk alik et all (120151 ) 
and lSzabo et ali (120141 4 for Kepler and CoRoT stars, respec¬ 
tively. We split the data sets into 1.5-day (EPIC 60018238 
and 60018678) and 1-day (EV Psc) bins and calculated the 
Fourier-parameters in the bins. As Figure [5] illustrates, all 
stars show some variation, even on this short time scale. It 
is of course possible that there are actually multiple unre¬ 
solved, separate peaks in the spectra but the considerably 
longer observations in the prime mission showed no signs of 
closely-spaced peaks in other RRc stars. 

We carried out the same analysis on the RR.d stars as 
well, but we did not detect significant amplitude variations 
in the fx mode. However, we found the uncertainties of the 
amplitude and phase data to be much higher than in the 
RRc stars, possibly because of the more complicated fre¬ 
quency content. Variation in the fx mode of t h e star CoRoT 
ID 0101368812 was detected bv ISzabo et al.l (120141 1. there¬ 
fore we expect that the science data of K2 will be sufficient 
to probe the variation of these modes in other RRc and RRd 
stars. 

3.3 A strongly modulated RRc star 

The fourth star, CSS J235742.1-015022 is located within the 
subframe of an engineering target star, EPIC 60042292, but 
it has no EPIC identification itself. The K2 observations re¬ 
vealed that it shows very strong amplitude and phase mod¬ 
ulation, as seen in Figure [6] The Blazhko effect exists in 
first-overtone RR Lyrae stars but i t is less commo n than 
in their fundamental-mode siblings. iMizers H (1200311 found 
a 10 per cent inciden c e rat e in the OGLE-II bulge data 
while (Nagy fe Kovacd (|2006i l derived 7 .5 per cent for the 
Large Magellanic Cloud based on MACHO data, although 
both numbers should be considered as lower limits only. Re¬ 
cent surveys, such as the Northern Sky Variability Survey 
and the All-Sky Aut omated Surv ey id entified several modu - 
lated field RRc stars dWils et al.ll2006l ; ISzczvgiel fc Fabrvckvl 


Figure 6 . Top: K2-E2 light curve of the modulated RRc star. 
Bottom: the Catalina Sky Survey (black) and K2 data (orange, 
shifted to the CSS brightness level), and folded with the modu¬ 
lation period. 

1200711 . The Blazhko Stars in the Gal actic Fi eld database 
currently lists 55 modulated RRc stars (jskar kall2013h . First- 
overtone stars with obvious Blazhko effect have not been o b¬ 
served before from space, though we note that ISzabo et alj 
(l2014ll found a very small side peak in one of the CoRoT 
RRc stars (ID 0105036241) that might be connected to am¬ 
plitude and/or phase variations. 


3.3.1 Modulation properties of the RRc star 

The K2-E2 data covers less than a modulation cycle, there¬ 
fore we analysed the CSS observations to determine the 
modulation period. We identified a single potential side 
peak, fi+fm, and determined a tentative modulation period 
of Pm = 17.07 d. Comparison with the K2-E2 data revealed 
that this is indeed the correct value and it has been stable 
for more than 8 years, over the duration between the start of 
CSS and the K2 observations (lower panel of Figure [6]). The 
detection of one modulated star out of eight RRc targets in 
the current Kepler and I<2 sample, agrees with the findings 
of previous studies and confirms the lower occurrence of the 
Blazhko effect among first-overtone stars. Adding the two 
CoRoT stars, the ratio of modulated stars falls somewhere 
between 10-20 per cent, depending on whether we consider 
CoRoT ID 0105036241 modulated or not. More observations 
will be required of course to determine the occurrence rate 
more precisely. K2 is perfectly suited to provide the answer: 
we expect it to observe approximately 50-100 galactic RRc 
stars over the course of the mission. 

We divided the data of CSS J23542.1-015022 into twelve 
0.75-day long segments and calculated the amplitudes and 
phases of the /1 frequency peak in each bin to investigate 
the modulation properties of the star. The results are shown 
in the lower panels of Figure [7] The pulsation amplitude 
changes rather symmetrically, but the descending branch 
of the phase variation is much steeper than the ascend¬ 
ing branch. In a phase-amplitude diagram, the modulation 
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CSS J235742.1 -015022 - K2-E2 data spectrum 
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these problems. Further observations can also provide clues 
to whether the modulation in RRab and RRc stars is the 
same mechanism or is caused by different physical processes. 


4 FUNDAMENTAL-MODE STARS 

We found 27 RRab-type stars in the K2-E2 data. Identi¬ 
fications and folded light curves are presented in Table [5] 
and Figure [5] The brightness range of the stars are similar 
to the original Kepler RR Lyrae sample (excluding RR Lyr 
itself), between Kp = 12.199 and 17.462 magnitudes. We 
were able to extract high-precision photometry for almost all 
stars. The image of EPIC 60018780 is affected by dark CCD 
columns (negative video crosstalk). We extracted the flux 
from the remaining pixels, but because the pointing drifts 
affect the dark columns and the star differently, small jumps 
can be seen in the data, especially in the first part of the 
observations. Given the limited resolution of Kepler , some 
stars are affected by blending. We see the strongest con¬ 
taminations from another nearby object in EPIC 60018773 
and 60018737. A more detailed evaluation and possibly im¬ 
age subtraction methods will be required to extract better 
photometry for these stars. 


Figure 7. Top: Fourier-spectra of the modulated RRc star, CSS 
J23542.1-015022 (within the mask of EPIC 60042292), original 
with thin grey line, prewithened with thick black line. Middle: 
residual spectra after standard (thin grey) and time-dependent 
(thick black) prewithening with the pulsation frequencies. The 
dashed arrows marks the expected position of the fx mode. Bot¬ 
tom: temporal variation of the amplitude and phase of the /i 
frequency, folded with the modulation period. 


progresses in the clockwise direction, opposite the direction 
that was observed in a nother modulated RRc star, TV Boo 
llSkarka fc Zeidall2013i ). 

3.3.2 Search for additional modes 

The Fourier spectrum of the star is shown in the up¬ 
per and middle panels of Figure d We found no signs of 
additional frequencies in the star. As the K2-E2 data is 
shorter than the modulation period, the modulation side- 
peaks cannot be resolved with conventional prewithening 
(constant Fourier parameters) properly. We applied an al¬ 
ternative, time-dependent prewithening method that sub- 
tracts non-stat i onary signals from the data, developed by 
iMoskalik et all J2015I ). While this method resulted in much 
lower residuals, we still could not identify additional modes 
down to 0.4 millimag amplitude. Interestingly, ground-based 
observations of two other modulated RRc stars did not re¬ 
veal any additional modes either, down to about 5 and 
2 millimagnitudes for LS Her and TV Boo, resp ectively 
llWils. Kleidis fe Broensll2008l : ISkarka fe Zeidall2013l ). 

The fx mode seem to be a common occurrence in nor¬ 
mal first-overtone stars that show only low-level tempora l 
variability but not the Blazhko effect (IMoskalik et al.|[2015l ). 
The potential lack of the fx mode in strongly modulated 
RRc stars may raise several questions about mode excita¬ 
tion and the physical origin of the modulation. The K2 mis¬ 
sion is ideally suited to expand our understanding about 


4.1 The Blazhko stars 

Although the 8.9-day long time base is rather short, we iden¬ 
tified amplitude variations in several stars. 11 out of the 27 
stars show significant residual power in their Fourier spectra 
after the removal of the pulsation frequency and its harmon¬ 
ics. While the data are too short to resolve the modulation 
side peaks, the light curves show smoothly changing pulsa¬ 
tion amplitudes. Therefore it is safe to assume that these 
stars are modulated. In 13 cases we found no signs of am¬ 
plitude and/or phase change either in the spectra or in the 
amplitude variation data. Analysis of 3 stars led to incon¬ 
clusive results. Therefore at least 40.7 per cent of the RRab 
stars seem to be modulated, but the number can be as high 
as ~ 50 per cent, in accordance with previous ground-base d 
and Kepler findings (ljurcsik et al.ll2009l ; lNemec et al.ll2013l ). 
This is an impressive result if we consider that Blazhko pe¬ 
riods are usually much longer than 9 days (in the Kepler 
sample they ranged from 27 days to two years). 

In three cases, data problems prevented us from con¬ 
firming or excluding the presence of modulation. EPIC 
60018737 and 60018773 suffer from blending. In both cases, 
the Fourier spectra did not reveal any clear and significant 
modulation sidelobes. In fact, slight temporal variation can 
be detected in the amplitude of the /1 peak, but the light 
curve itself is distorted and includes additional variation too. 
EPIC 60018780 is severely impacted from the negative video 
crosstalk which resulted in jumps in the light curve and el¬ 
evated noise level in its Fourier spectrum. 

f.1.1 Modulation properties of RRab stars 

The K2-E2 data in themselves are too short to determine 
the modulation periods and amplitudes of the stars. While 
some of the brighter non-modulated stars were covered by 
multiple ground-based surveys, for the Blazhko stars we had 
to rely on the CSS data alone. The Catalina data has very 
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Figure 8. Folded light curves of the 27 RRab stars in the K2-E2 sample, based on the periods presented in Table [5] In black, roman 
labels: non-modulated stars; in orange, and marked with asterisks: Blazhko stars; in grey, italic labels: uncertain. 


sparse sampling in this area with 2-300 data points scattered 
over 7-8 years. Each field is covered 3-4 times per night and 
fields repeat after one or two weeks. This kind of sampling is 
well suited for asteroid hunting but less useful to determine 
the Blazhko periods for RR Lyrae stars. The precision of 
individual CSS data points is 0.05 mag for bright targets 


(~ 13 mag) but reaches 0.1 magnitudes for targets fainter 
than 17 mag that is well below the (sub)mmag accuracy of 
the K2 data. 

We searched for possible modulation sidepeaks in the 
Fourier spectra of the CSS data. In some cases the beating 
between the pulsation period and the sampling resulted in 
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Table 5. Fundamental-mode RR Lyrae stars observed during the Two-Wheel Concept Engineering Test run of the K2 mission. Photo¬ 
metric [Fe/H] indices are included as well. The average uncertainty of the indices is ±0.1 dex. 


EPIC ID 

RA 

(deg) 

Dec 

(deg) 

Kp 

(mag) 

Period 

(d) 

Name 

[Fe/H] 

CSS ID No. 

Pm 

(d) 

Addtl. 

modes 




(a) Non-modulated RRab stars 





60018227 

0.83750 

±3.39833 

13.720 

0.579041 

ASAS J000321±0323.9 

-1.30 

- 

- 

- 

60018644 

354.15513 

-2.21173 

12.199 

0.645043 

ASAS J233637—0212.7 

-1.88 

1001126012320 

- 

g-mode? 

60018650 

356.29939 

-5.20608 

13.636 

0.630289 

NSVS 11906749 

-1.95 

1004127008648 

- 

- 

60018663 

355.36115 

-5.54157 

14.221 

0.468800 

NSVS 11904882 

-1.34 

1004127002457 

- 

- 

60018669 

351.61480 

-2.97015 

14.490 

0.465721 

CSS J232627.6—025814 

-1.25 

1004126052487 

- 

- 

60018670 

357.53961 

±1.02409 

14.501 

0.669217 

NSV 26151 

-1.29 

1101128018965 

- 

- 

60018679 

1.1001 

-6.85657 

14.888 

0.494723 

CSS J000424.0—065123 

-1.30 

1007001029423 

- 

- 

60018708 

356.48300 

-7.47260 

15.813 

0.590650 

CSS J234555.9—072821 

-1.63 

1007126018552 

- 

- 

60018722 

359.57559 

-8.71836 

16.133 

0.482625 

CSS J235818.1-084306 

-1.29 

1009127050514 

- 

- 

60018732 

352.76588 

±1.06081 

16.292 

0.498292 

CSS J233103.7±010339 

-1.56 

1101126021705 

- 

- 

60018752 

352.46439 

±1.65528 

16.710 

0.594034 

CSS J232951.4±013919 

-1.26 

1101126033293 

- 

(PD?) 

60018755 

1.94723 

±3.26733 

16.740 

0.681966 

CSS J000747.3±031602 

-1.78 

1104001007772 

- 

- 

60018779 

351.77498 

-2.31793 

17.415 

0.465741 

CSS J232705.9-021904 

-1.33 

1001126009779 

- 

- 




(b) Modulated RRab stars 





60018656 

355.98041 

-4.05090 

13.956 

0.565113 

CSS J234355.2—040303 

-1.54 

1004127030869 

- 

- 

60018657 

353.86545 

-5.72775 

14.054 

0.543388 

CSS J233527.6—054340 

-1.38 

1007125047274 

84.3 

PD, F2 

60018673 

357.75112 

±4.87943 

14.561 

0.620309 

CSS J235100.2±045245 

-1.49 

1104128039542 

- 

F2 

60018676 

353.11519 

-5.39302 

14.633 

0.699815 

CSS J233227.6—052334 

-1.68 

1004126005022 

65.5 

- 

60018709 

352.8671 

-2.62876 

15.815 

0.563396 

CSS J233128.1-023743 

-1.19 

1001126004067 

- 

- 

60018742 

1.90229 

-8.24558 

16.460 

0.594456 

CSS J000736.5—081444 

-1.50 

1007001004138 

- 

- 

60018743 

352.99328 

-4.92403 

16.481 

0.493351 

CSS J233158.3-045526 

-1.75 

1004126013659 

- 

PD 

60018750 

356.47678 

-7.05360 

16.658 

0.548491 

CSS J234554.4—070312 

-1.36 

1007126026371 

- 

- 

60018763 

357.49761 

±4.62398 

16.900 

0.586196 

CSS J234959.4±043726 

-1.19 

1104128034480 

- 

- 

60018771 

357.55888 

±4.24208 

17.065 

0.553989 

CSS J235014.1±041431 

-1.22 

1104128027645 

- 

- 

60018778 

355.94061 

±1.86808 

17.376 

0.620342 

CSS J234345.7±015205 

-2.00 

1101127037513 

- 

- 




(c) Possibly modulated RRab stars 





60018737 

355.31566 

-4.21736 

16.376 

0.635342 

CSS J234115.7-041302 

-1.26 

1004127027477 

- 

- 

60018773 

359.53266 

±3.27389 

17.094 

0.595290 

CSS J235807.8±031626 

-1.30 

1104128008817 

- 

- 

60018780 

1.56091 

-4.50438 

17.462 

0.505453 

CSS J000614.6—043015 

-1.08 

1004001019454 

- 

- 


modulation-like variations. To exclude them, we compared 
the modulation phases between the CSS and K2 light curves. 
We were able to determine the Blazhko periods for two stars, 
EPIC 60018657 (Pm = 84.3 ± 0.3 d) and 60018676 (Pm = 
65.5±0.3 d). Despite the poor coverage, the final light curves 
in Figure O line up quite well. Interestingly, this method 
failed for two stars that show strong amplitude variations, 
EPIC 60018742 and 60018743. We suspect that these stars 
have shorter modulation periods than the two we presented 
above and the CSS observations are simply too sparse in 
time to identify those periods properly. 


Finally, EPIC 60018709 shows variations in its light 
curve shape and pulsation amplitude. After prewithening 
with the pulsation frequency and its harmonics, however, 
very little evidence can be found for side peaks anywhere 
but the 3/i frequency value. Examination of the CSS data 
did not reveal any significant modulation side peaks, but 
the folded light curve indicates some amplitude variation 
around maximum light. Therefore we classihed this star as 
a modulated RRab variable. 


EPIC 60018657 - CSS J233527.6-054340 
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Figure 9. Folded Catalina Sky Survey and K2-E2 light curves of 
two modulated RRab stars where we were able to estimate the 
Blazhko periods. 
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EPIC 60018657 - prewithened spectrum 



EPIC 60018644 - prewithened spectrum 



Figure 10. Top: residual Fourier spectrum of the modulated star 
EPIC 60018657 after prewithening with the pulsation frequency 
and its harmonics (grey dashed lines). Identified or suspected ad¬ 
ditional peaks are labelled. In this case the Blazhko period is 
long (84.3 d), so the star shows little amplitude and phase vari¬ 
ation over the K2 observations leading to very small modulation 
sidepeaks in the spectrum. Bottom: residual spectrum of EPIC 
60018644, revealing the peaks corresponding to the g-mode sus¬ 
pected in the star. Thin solid and dashed arrows indicate the com¬ 
bination frequencies, /o + nf g and n/o — f g , respectively. Some 
labels are not shown in the figure to avoid clutter. 

4-1.2 Additional modes in RRab stars 

One of the major discoveries of Kepler, regarding 
fundamental-mode RR Lyrae star s, was the identification o f 
low-amplitude additional modes (Benko et al. 20ld . 120141) . 
With the help of model calculations, we were able to trace 
back the origin of period doubling to one of these modes , 
the ninth overtone (ISzabo et al.l 1 20 id : iKollath et al.l 1201 ill . 
We note in passing that although the ninth overtone itself 
has very low amplitude and is undetectable in the K2-E2 
sample, period doubling itself can be still observed since it is 
not a classical beating between modes but a purely dynami¬ 
cal phenomenon. The strong resonance itself destabilizes the 
limit cycle of the fundamental mode that bifurcates into a 
new, period-doubled limit cycle that takes two pulsation pe¬ 
riods to repeat. This in turn creates a series of half-integer 
peaks in the Fourier spectrum but those do not correspond 
to actual pulsation modes themselves. 

We searched thoroughly for additional modes in the 
K2-E2 data set. For classical variables, the widely used cri¬ 
terion of accepting a peak is S/N > 4. In the long data 
sets available from C 0 R 0 T and the Kepler prime mission, 
however, combinations with fo were always observed for the 
detected additional peaks. Therefore, we also considered a 
peak (or a series of peaks) reliable and of stellar origin when 
it formed combination peaks with the fundamental mode, 
e.g., a regular pattern was detectable in the spectrum even 
if the combinations had a S/N ratio between 2 and 4. Four 
authors searched for periods independently and the results 
were compared to each other: we only accepted peaks that 
were detected by at least two authors. The results are in¬ 
cluded in Tableland we detail them below. 
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Figure 11. Fourier spectrum of the RRab star EPIC 60018673 
(top), and the residual spectra. Middle: residual of the standard 
prewithening with the main frequency and its harmonics: the 
strong, unresolved residual power hinders the detection of the 
additional modes. Bottom: residual of the time-dependent pre¬ 
withening. The /2 peak is unambiguously detected. 
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Figure 12. Fourier spectrum of the RRab star EPIC 60018752 
(top), and the residual spectrum (bottom). A single half-integer 
peak (1.5/o) is tentatively detected, suggesting that period dou¬ 
bling might occur in the pulsation. 


We identified the same instrumental frequency peak 
that we found in EPIC 60018238, / corr ~ 8.16 d _1 , in two 
RRab stars, EPIC 60018227 and 60018644, and excluded it 
from the analysis. 

Despite our efforts, we found period doubling in only 
two stars, EPIC 60018657 and 60018743, and we have ten¬ 
tative detection in one more faint target (see Figure fl2l) . We 
now know, however, that the amplitude of period doubling 
is very variable, often on short time scales. Long-term obser¬ 
vations provided by C 0 R 0 T and Kepler can help the detec¬ 
tion by accumulating multiple high-amplitud e phase s and 
lowering the noise in the Fourier-spectra feenko et al.ll2014l : 
ISzabo et al.1120141 ). On the other hand, short-term measure¬ 
ments, such as the K2-E2 can miss these signals entirely. 

The Fourier spectrum of EPIC 60018657 revealed not 
only half-integer peaks ((2n + l)/2 fo), the tell-tale signs of 
period doubling, but other additional peaks as well. As the 
residual spectrum in Figure flOl illustrates, the /2 = 3.14198 
d -1 mode and its linear combinations with the n/o frequen- 
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cies can be easily detected. Its frequency ratio, / 0//2 = 
0.5852, is very similar to the ratios observ ed in other Blazhko 
stars in the Keple r and C 0 R 0 T samples l|Benko et, al.ll2014l : 
ISzabo et al.ll20l3 l. This mode is hypothesized to be the sec¬ 
ond radial overtone, although it was not reproduced by non¬ 
linear hydrodynamic models so far. The reason of its absence 
could be that it is in fact a non-radial mode close to the fre¬ 
quency of the second overtone. It can also be hypothesized 
that the excited mode exists in a parameter space that ID 
models have avoided so far, similar to the discovery of pe¬ 
riod doubling and the l ow-am plitude first overtone in some 
RRab models (IMolnar et al.ll2012r ). We note, however, that 
the mode was not detected during the extensive survey of 
those models. Some excess amplitude still remains in the 
spectrum after the removal of the half-integer and /2 peaks, 
but the dataset is too short to resolve it into further addi¬ 
tional modes. We detected the peak firmly in one other 
star, EPIC 60018673, at /2 = 2.7380 d -1 , with a period ratio 
of / 0//2 = 0.5881. This star shows stronger side peaks but 
we were able to detect /2 with the time-dependent Fourier 
method, despite its small amplitude (1.1 mmag, Figure fill) . 

Additionally, we identified an intriguing periodicity in 
EPIC 60018644. Although the peak at / = 1.0774 d -1 itself 
has a signal-to-noise ratio of 3, we identified multiple signif¬ 
icant peaks related to it (harmonics and combinations with 
/o), therefore we accept this as a potential detection of a sep¬ 
arate mode in the star. What makes this mode particularly 
interesting is that the period is longer than the fundamental 
mode itself, with a P/Po = 0.695 ratio, so it must be a g- 
mode that are always non-radial. Similar modes have been 
identified recently in first-overtone stars observed by Kepler 
and OGLE. The additional mode found in multiple OGLE 
stars has a period ratio of ~ 0.92 — 0.95 (scaled to the ex¬ 
pected period of the fundamental mode) that is very differ¬ 
ent from our value (iNetzel. Smolec fc Dziembowskil l2015lh 
On the other hand, Kepler RRc stars show a multitude 
of modes between sc aled period ratios from 0.38 to 0.97 
llMoskalik et al.ll2QI5l ). 

As the last column of Table 0 indicates, we did 
not detect unambiguous additional p-modes in the non- 
modulated RRab stars, in accordance with previous results 
dBenko et al.ll2C)i~4l : ISzab6 et al.ll2014l l. although a potential 
g-mode was identified in one star. Given the shortness of the 
data sets, however, both the tentative mode detections and 
the non-modulated status of the stars can be erroneous. 
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Figure 13. Fundamental periods versus <j >31 Fourier phase rela¬ 
tions of the RRab stars. The relation of the two measures provide 
the photometric metallicities of the fundamental-mode stars that 
are included in Table [5] Modulated stars are indicated with addi¬ 
tional boxes around the squares. Dashed li nes corresp ond to va r- 
ious iso-metallicity lines, based on Eq. 3 of lNemec et ali (2013 '). 


The accuracy is likely more limited for the modu¬ 
lated stars where the phase relation valu es change over 
the Blazhko cycle, but it is still acceptable. ISmoled i 20lV» i 
showed that filtering out the modu lation side-pea ks can lead 
to reasonably good [Fe/H] values. iNemec et alj (l2013li also 
investigated the relation between the spectroscopic and pho¬ 
tometric [Fe/H] determinations and found that the values 
are consistent for most Blazhko stars. Only three extremely 
modulated stars, less than 10 per cent of the overall sample, 
showed significant discrepancies. Based on the original Ke¬ 
pler sample, we can estimate that about 90 per cent of the 
K2-E2 stars should closely follow the relation we used. 

We consider the limited modulation phase coverage the 
strongest source of uncertainty in our study. Both the pul¬ 
sation period and cf >31 change duri ng the modul a tion c ycle: 
based on the values presented by INemec et alj (|2013l l , we 
estimate the overall uncertainty of the [Fe/H] values for the 
Blazhko stars to be ±0.15 — 0.2 dex. 


4.2 Photometric [Fe/H] indices 

Light curve parameters, especially the (j >31 Fourier phase 
relat ions of RRab stars (</31 = (f >3 — 3 </>i, as defined by 
ISimon fc Led ll98ll l ca n be used to calcu late photomet¬ 
ric metallicity indices. INemec et al] d2013il compared the 
Fourier parameters of the original Kepler RR Lyrae sample 
with indices derived from spectra obtained with the CFHT 
and Keck te lescopes. W e applied the same relation (Eq. 3 
from INemec et aT1l2013l l to calculate the photometric [Fe/H] 
indices of the K2-E2 stars. The results are plotted in Figure 
m and included in Table [5] The K2-E2 sample is moder¬ 
ately metal-poor, with [Fe/H] valu es rang ing fro m -1.08 to 
2.00. Based on the experiences of lNemec et al.1 d2013l ~). we 
estimate the accuracy of the indices to be ± 0.1 dex for non- 
modulated stars. 


5 CONCLUSIONS AND FUTURE PROSPECTS 

The Kepler space telescope observed 1952 targets in Pisces 
for 8.9 days during the K2 Two-Wheel Concept Engineer¬ 
ing Test run. We identified 32 RR Lyrae stars among the 
targets and one more as a background object in the target 
pixel mask of another star. We extracted the light curves 
from the target pixe l files with the PyKE software tool 
([Still fe Barclavll2012lh The analysis of the photometric data 
yielded the following results: 

• We identified 2 double-mode (RRd), 4 first-overtone 
(RRc) and 27 fundamental-mode (RRab) stars in the held. 
Mos t of the stars were discovered by the Catalina Sky Sur¬ 
vey dDrake et al.ll20ljl . 

• The fx additional mode was detected in the two RRd 
and the three non-modulated RRc stars, although with vary- 
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iiig degrees of significance. The mode appears to be period- 
doubled and varies over time in three out of five stars. 

• One of the RRc stars, CSS J125742.1-015022, turned 
out to be a modulated first-overtone star. This is the first 
time that a Blazhko-RRc star was observed by a space 
photometric mission. The modulation period of the star is 
17.07 d and both the pulsation amplitude and phase show 
strong variations. It is the only RRc star where no additional 
modes were found. 

• Indications of a potential g -mode was found in one 
RRab star, EPIC 60018644, at a period ratio of Po/P 9 = 
0.695. 

• Even from such short dataset, we found indications that 
at least 11 out of 27 RRab stars (41 per cent) may be mod¬ 
ulated. We detected period doubling and the (presumed) 
second overtone in multiple stars. Based on the Fourier- 
parameter relations, the [Fe/H] values of the stars fall be¬ 
tween — 1.08 and —2.00 dex. 

• As a by-product of the target search for the K2 cam¬ 
paigns, we re-examined the ‘RRd’ and ‘Blazhko’ classifica¬ 
tions of the Catalina Periodic Variable Star Catalog and 
identified 165 bona-fide double-mode stars throughout the 
sky, 130 of which were previously unknown. For a brief de¬ 
scription, see Appendix 0 

The results promise a bright future for the K2 mission. 
Even from these very short datasets, we were able to identify 
all new phenomena that space-based photometric missions 
discovered in RR Lyrae stars so far. The 75-day long cam¬ 
paign data will provide much better frequency resolution 
and refined light curve extraction methods will improve the 
photometric precision. 

Losing the ability of Kepler to observe stars for several 
years was unfortunate and will be very hard to reproduce. 
The silver lining of K2 mission, however, is its flexibility to 
observe many different stars. This single field contained 33 
RR Lyrae variables. Based on our estimates and target se¬ 
lection so far, this number of expected RR Lyrae stars will 
vary b etween approximatel y 10 and 200 for the forthcoming 
fields jMoInar. Plachv fe Szaboll2014l L With a sustained ob¬ 
serving program, we will be able to answer new questions the 
previous missions could not: what is the true ratio of modu¬ 
lated RRc stars? Are all non-modulated RRab stars strictly 
periodic? Is the fx mode present in all RRd and RRc stars? 
Finally, K2 will be able to detect RR Lyrae stars in glob¬ 
ular clusters, halo streams, and hopefully in nearby dwarf 
galaxies as well, broadening the scope of the mission even 
further. 
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(l2014ll . The stars we found are distributed between declina¬ 
tions —30° and 60°, except for a ~ 30-degree band along 
the Galactic plane (Figure TAlll . 

Double-mode RR Lyr stars in t he Catalina sampl e 
were separated from other RR Lyr by iDrake et alj (120141 s ), 
based on the visual examination of the light curves. No 
prewhitening was performed, thus, such classification can 
be uncertain as indicated by the misclassification of the two 
stars in the K2-E2 field. As the rare RRd stars are valu¬ 
able targets, we searched for bo na-fide double- mode pul- 
sators among the stars that IDrake et alj (120141 s ) classified 
as ‘RRd’ or ‘Blazhko’ stars. The photometri c data were 
downloaded from http://catalinadata.org/ (I Drake et al.1 
120091 s ) . We found 461 stars classified as RRd and 70 stars 
classified as Blazhko for which at least 100 data points 
were obtained. Double mode pulsations were found in 161 
and 5 stars, respectively. This means that only 35 per cent 
of the stars classified as ‘RRd’ in Catalina catalogue are 
in fact double-mode pulsators. We checked that the frac¬ 
tion of confirmed RRd stars does not change if we limit 
to stars with at least 300 data points. Most of the stars 
that were erroneously identified as ‘RRd’ in the CSS vari¬ 
able catalogue actually show Blazhko-type modulation. One 
star turned out to have two different CSS identifications 


(CSS_J104003.0+414504 = CSS.J104003.2+414503), bring¬ 
ing the total number of confirmed RRd stars to 165. We 
checked the literature for these stars and found that 35 of 
them were previou sly identified, most in the LINEAR sur¬ 
vey data (29 stars, Poleski 2014). Six stars were previousl y 
found by [ciementiniet^! ( 2000()jjBernaTd& ii Wn^ 


— - —- -— - 1 1 • mii .ti i - * • 7 

Szczy giel fc Fabrvckv (2003); Hambsch fe Wild ( 200Slb and 
Wild 1 20101 s ) . The number of confirmed RRd stars is increased 
by 130. Table IaTI presents all confirmed RRd stars and Ta¬ 
ble [A2] lists their pulsation parameters. Figure [2] compares 
the positions of newly identified stars on Petersen diagram 
with RRd stars from various other sources. 

We also checked CSS photometry for the two K2- 
E2 RRd stars. After removing outlying data points from 
Catalina data we were able to confirm double-mode pulsa¬ 
tions in both stars. This illustrates that more double-mode 
stars are hidden in Catalina lists of the RRab and RRc stars. 
However, those lists are much longer, with tens of thousands 
and more than 5000 stars, respectively, and will require fur¬ 
ther work to classify. 


APPENDIX A: DOUBLE-MODE STARS IN THE 
CATALINA SKY SURVEY 

As a side-product to the K2 investigations we also identi¬ 
fied 130 new RRd stars throughout the sky. These are very 
valuable targets, because their physical parameters can be 
estimated from hydrodynamic model calculations, as sec¬ 
tion [3TT7T] illustrates. The Catalina Sky Survey provided us 
perhaps with the richest database of RR Lyrae stars for tar¬ 
get selection for all K2 fields that are not too close to the 
Galactic plane. As usually only 1-5 RRd stars can be ob¬ 
served in a K2 campaign, proper identification is especially 
important. This work is also the continuation of the search 
for bona-fide RRd stars in the LINEAR sample bv IPoleskH 
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Table Al. Sample table of the confirmed double-mode RR Lyae stars in the Catalina Surveys Periodic Variable Stars Catalog. The 
list includes 165 stars, based on the examination of the light curves of the stars originally labelled as ‘RRd’ or ‘Blazhko’ in the Catalog, 
along with the two K2 stars. Full table is available in the online version. 


CSS name 

RA 

(deg) 

Dec 

(deg) 

Brightness 

(mag) 

CSS ID 

Other ID 

CSS.J001420.8+031214 

3.586833 

+3.20389 

17.45 

1104002007409 

— 

CSS.J001724.9+200542 

4.353750 

+20.09506 

16.64 

1121002007726 

— 

CSS_J001812.9+210201 

4.554042 

+21.03375 

14.54 

1121002027610 

— 

CSS-J001836.2-191522 

4.651083 

-19.25619 

16.84 

1018002009928 

— 

CSS_J003359.4+022609 

8.497833 

+2.43583 

15.87 

1101004049971 

— 

CSS.J004804.2+365310 

12.017833 

+36.88625 

16.60 

1138004006108 

— 

CSS_J011032.5+215708 

17.635500 

+21.95233 

16.13 

1121006050243 

— 

CSS-J011824.3-123318 

19.601583 

-12.55525 

17.72 

1012007028731 

— 

CSS-J014305.3+010549 

25.772167 

+1.09708 

17.04 

1101010021310 

— 

CSS_J082622.9+282405 

126.595625 

+28.40142 

16.84 

1129040004714 

LINEAR 658512 


Table A2. Sample table of the pulsation parameters of the confirmed double-mode RR Lyae stars. FM: fundamental mode, Ol: first 
overtone. Full table is available in the online version. 


CSS name 

FM period 
(day) 

FM period 
uncert. 

FM ampl. 
(mag) 

Ol period 
(day) 

Ol period 
uncert. 

Ol ampl. 

(mag) 

Per. ratio 

(Pi/Po) 

CSS-J001420.8+031214 

0.51930994 

0.00000166 

0.23 

0.38711458 

0.00000056 

0.39 

0.74544034 

CSS-J001724.9+200542 

0.48002554 

0.00000174 

0.29 

0.35712249 

0.00000074 

0.41 

0.74396560 

CSS.J001812.9+210201 

0.55837503 

0.00000198 

0.19 

0.41615750 

0.00000064 

0.34 

0.74530106 

CSS-J001836.2-191522 

0.46299846 

0.00000198 

0.28 

0.34403253 

0.00000067 

0.40 

0.74305329 

CSS.J003359.4+022609 

0.48368923 

0.00000208 

0.21 

0.36009918 

0.00000087 

0.29 

0.74448459 

CSS-J004804.2+365310 

0.46874750 

0.00000070 

0.38 

0.34831989 

0.00000079 

0.37 

0.74308640 

CSS.J011032.5+215708 

0.47302180 

0.00000138 

0.33 

0.35178257 

0.00000051 

0.40 

0.74369209 

CSS-J011824.3-123318 

0.52620846 

0.00000320 

0.27 

0.39247062 

0.00000147 

0.38 

0.74584628 

CSS-J014305.3+010549 

0.47555267 

0.00000161 

0.32 

0.35372017 

0.00000058 

0.40 

0.74380861 

CSS.J082622.9+282405 

0.47390930 

0.00000138 

0.22 

0.35260694 

0.00000044 

0.35 

0.74403887 


CSS RRd stars 


known • 



Figure Al. Distribution of the CSS double-mode RR Lyrae stars 
in the sky, in Equatorial coordinates. Blue circles are the new 
discoveries, orange dots are the already known stars. 















